Random deposition of particles of different sizes 
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We study the surface growth generated by the random deposition of particles of different sizes. 
A model is proposed where the particles are aggregated on an initially flat surface, giving rise to a 
rough interface and a porous bulk. By using Monte Carlo simulations, a surface has grown by adding 
particles of different sizes, as well as identical particles on the substrate in (1 + 1) dimensions. In 
the case of deposition of particles of different sizes, they are selected from a Poisson distribution, 
where the particle's sizes may vary by one order of magnitude. For the deposition of identical 
particles, only particles which are larger than one lattice parameter of the substrate are considered. 
We calculate the usual scaling exponents: the roughness, growth and dynamic exponents a, j3 and 
z, respectively, as well as, the porosity in the bulk, determining the porosity as a function of the 
particle size. The results of our simulations show that the roughness evolves in time following three 
different behaviors. The roughness in the initial times behaves as in the random deposition model. 
At intermediate times, the surface roughness grows slowly and finally, at long times, it enters into the 
saturation regime. The bulk formed by depositing large particles reveals a porosity that increases 
very fast at the initial times, and also reaches a saturation value. Excepting the case where particles 
have the size of one lattice spacing, we always find that the surface roughness and porosity reach 
limiting values at long times. Surprisingly, we find that the scaling exponents are the same as those 
predicted by the Villain-Lai-Das Sarma equation. 

PACS numbers: 68.43.Jk ; 68.35. Ct ; 81.15.Aa ; 02.50.-r 



I. INTRODUCTION 

The research on the geometric properties of growing 
surfaces is one of the most important in the field of the 
non-equilibrium statistical physics, not only because it is 
a challenge for the theoretical physicist to model these 
properties, but especially for the intrinsic experimental 
interest in shaping the surfaces with desired purposes 
0, 0. The random deposition is the simplest known 
model, where particles are aggregated onto an initially 
flat substrate. Because lateral correlations among the 
deposited particles are completely neglected, the contin- 
uous and discrete atomistic versions of the model have 
exact solutions. If a surface relaxation mechanism is 
allowed to the deposited particles in the random de- 
position model, height-height correlations naturally ap- 
pear. Although the corresponding discrete model did not 
present an exact solution, it can be described by the lin- 
ear Edwards- Wilkinson equation (EW) @, which is ex- 
actly soluble. A generalization of the (EW) equation was 
proposed by Kardar, Parisi and Zhang in 1986, including 
a non-linear term that accounts for the lateral growth 
of the interface. Although the growth process is a local 
phenomenon, the lateral growth is related to the spread- 
ing of the heigth fluctuations along the surface. This is 
characterized by a correlation length that increases with 



time, and it reaches a maximum value corresponding to 
the linear dimension of the substrate. A perpendicular 
correlation length is also defined that is related to the 
fluctuations in heigth along the growth direction. The 
non-linear KPZ equation [4j is very useful to describe 
porous deposits, as those ones generated by the ballistic 
deposition of particles @. 

The morphology of the surface is described by the in- 
terface width w(L,t), which characterizes the roughness 
of the interface, is defined by the rms deviation of the 
height h around is average value h, that is, 



w(L,t) = \ < (h~h) 2 >, 



(1) 



where L is the linear dimension of the substrate, t is the 
time elapsed after the growth start, the overbar indicates 
spatial average and the angular brackets mean configura- 
tional averages. It is well established that a large class of 
growth models follow the Family- Vicsek scaling relation 



w (L,t)~L a f( — ), 



(2) 
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where the scaling function f(x) is a constant when x is 
very large, and f{x) ~ x@ when x « 1. The exponent a 
characterizes the interface width, z is the dynamic expo- 
nent, while /3 is the growth exponent. These exponents 
are not independent and are related by a = j3z. 

Most of the studies performed in the area of surface 
growth focus their attention in the determination of these 
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exponents, and whether they fit to the Family- Vicsek 
ansatz. These calculations have been done analytically, 
by solving stochastic differential equations or employing 
mean field approximations, and through extensive use 
of Monte Carlo simulations 

[7 13]. In this work we use 
Monte Carlo simulations to study the surface growth in 
(1+1) dimensions due to the deposition of particles of dif- 
ferent sizes. There are in literature some studies where 
two or more different deposition models are combined 
0, [l5[ or two species of particles are deposited [l6l - fl8| 
in order to describe the time evolution of roughness in 
real systems. Here, we choose our particles to be de- 
posited from a modified Poisson distribution, with aver- 
age size equals to five lattice spacings and maximum size 
equals to nine. This type of distribution appears to be 
relevant in some ash particles deposition on the heat ex- 
change surfaces [IH, [2(| ■ We also consider in this work, 
the deposition of identical particles, which size is larger 
than one lattice spacing of the substrate. The deposition 
of particles proceeds as in the pure random deposition 
model, however, correlations between columns naturally 
appear due to the deposition of particles larger than one 
lattice spacing. The random deposition of particles of 
unit size generates a compact bulk and an infinitely large 
interface width. This means that for a linear substrate of 
size L, w(L,t) becomes infinitely large as the deposition 
time t goes to infinite. On the other hand, the deposition 
of particles larger than one unit produces a porous bulk. 
Then, we can determine the evolution of the porosity of 
the bulk simultaneously with the scaling properties of the 
surface. We say that the interface width satisfying Eq. 
(2) is self-affine, which means that rescaling part of the 
interface anisotropically we obtain an interface that is 
statistically indistinguishable from the whole one. In or- 
der to characterize the morphology of a rough interface, 
it is sufficient to know the value of the exponent a 
We show that the surface growth presents three different 
behaviors as a function of time. At the initial times, it 
behaves as in the usual random deposition of particles 
of unit size. Then, at intermediate times, when the lat- 
eral correlations between columns develop, it grows more 
slowly, and finally reaches a saturation regime at long 
times. Porosity attains the saturation regime faster than 
the interface width, and it presents large values even for 
the smallest particle of size two lattice spacings. Finite- 
size effects are not observed for this property, and for 
particles of size corresponding to four lattice units of the 
substrate, the porosity is already 85% of the maximum 
possible porosity, that is, P = 0.5. 

The paper is organized as follows. In Sec. II we present 
the model, the deposition rules and some details concern- 
ing the Monte Carlo simulations. In Sec. Ill, we present 
the results for the interface width as a function of par- 
ticles' sizes. Sec. IV contains the results for the scaling 
exponents when the deposited particles are identical, and 
Sec. V is devoted to the study of bulk porosity. Finally, 
in Sec. VI, we present our main conclusions. 



II. MODEL FOR DEPOSITION OF LARGE 
PARTICLES 

We propose a model where the particles are dropped 
randomly over a finite linear substrate, which is divided 
into cells of unit size. All the particles also have unit 
height and the columns where they land are increased 
by one unit. Particles of different sizes are selected from 
a modified Poisson distribution, with sizes changing by 
one order of magnitude. The most important difference 
between this model and the simple random deposition, 
where only particles of unit size are deposited, is that 
the present model naturally allows for correlations among 
the columns, that is, it leads to a lateral growth of the 
interface. We studied the effect of particles's sizes on the 
scaling exponents, as well as on the global porosity. In 
this study, we performed Monte Carlo simulations in two 
quite different conditions: I) adding particles of different 
sizes, and 2) adding identical particles. For the deposi- 
tion of identical particles, only particles which are larger 
than one lattice spacing of the substrate are considered. 
In the case of particles of different sizes, they are selected 
from a modified Poisson distribution. 

The Poisson distribution was modified in this work in 
order to take into account the discrete nature of the par- 
ticles and a maximum size of the particle to be deposited 
during our Monte Carlo simulations. The probability of 
a particle with size equal to n to be selected for a depo- 
sition trial is 



where n ranges from one to nine lattice spacings, which is 
the maximum permitted size for deposition in the present 
model. As we fixed x at the value five, the normalization 
factor a is equal to 0.96143. Therefore, this is a type of a 
modified Poisson distribution with a finite set of possible 
results. For a true Poisson distribution, where all values 
of n are possible, x should be the average value of the 
distribution and the normalization factor a = 1. This 
modified Poisson distribution is believed to be realistic to 
describe some ash deposits, where particles's sizes may 
vary by two or more orders of magnitude, and are well 
represented by a Poisson distribution [19j . 

When we perform simulations with particles of differ- 
ent sizes, the particles to be deposited are flat, height 
one, and may change from (1 X 1) to (9 X 1) in units of 
the lattice parameter of the linear substrate, according 
to the modified Poisson distribution. In the case of de- 
position of identical particles, we choose just one species 
of particle, with a defined size. The simulations in this 
work are performed in (1+1) dimensions, the resulting 
deposit is two-dimensional, and we always use periodic 
boundary conditions along the linear substrate. 

To add a particle on the substrate, a cell on the lattice 
is randomly selected. At this selected position we place 
the midpoint of the particle, trying to extend horizon- 
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FIG. 1: Illustration of the growth rules, where geometric re- 
strictions lead to the reflection of a particle by the surface. 
Particle "A" will be rejected, while particle "B" will be added 
to the surface. 



FIG. 2: Log-log plot of the roughness versus time for deposi- 
tion with particles ' size taken from the modified Poisson dis- 
tribution, Eq. (3). From bottom to top the lattice sizes are L 
= 1024, 2048, 4096, 8192. 



tally, at both sides of the initial position, this site to the 
full size of the particle. In the case of particles with even 
size, as 2, 4, 6 and 8 lattice spacings, we choose as their 
midpoint, the position 1, 2, 3 and 4, respectively. 

We allow the reflection of a particle by the surface and, 
in this case, the trial of deposition is lost, the particle is 
not incorporated to the substrate. The reflection is due 
to the geometric constraint properties of the surface. For 
each particle deposition attempt, we look at the number 
of vacant sites around a randomly selected position on 
the substrate. The number of vacant sites around this 
site determines the possible deposition or reflection of 
the selected particle. If the vacant region is smaller than 
the size of the selected particle, the particle will be re- 
jected. On the other hand, if the vacant region is larger 
or equal to the size of the particle, there are no geometric 
restrictions and the incident particle permanently sticks 
to that region. 

The growth rules are ilustrated in Fig. [1] The particle 
labeled by "A" will not be incorporated to the surface, 
while the particle labeled by "B" is added to the deposit. 
Under the "B" particle, a shadow zone is created, where 
no other particle can be aggregated. This is the mech- 
anism that leads to the formation of a porous structure 
inside the bulk. In order to determine the global porosity, 
we calculate the quantity 



P= , (4) 

Vpore ~\~ Vsolid 

where V pore is the volume occupied by the vacant sites, 
and V so ud is the proper volume occupied by particles. 
The vacant sites here mean only the empty sites not be- 
longing to the surface. The overhangs from the topmost 
layer are not considered vacant sites. In this work, we 
show that the deposition of particles larger than one lat- 
tice parameter of the substrate, generates a highly two- 
dimensional porous structure. The time evolution of the 
porosity and its dependence on the size of the deposited 
particles are presented in Sec. V. 



III. PARTICLES OF DIFFERENT SIZES 

In this Section, we focus our simulations on the esti- 
mate of the scaling exponents and determination of the 
universality class of the model, when the deposited parti- 
cles are selected from the probability distribution, Eq.(3). 
The usual scaling exponents: the growth, roughness and 
dynamic exponents, and also the global porosity in the 
bulk, are calculated for deposition of particles larger than 
one lattice spacing. Some earlier studies 12, 16], clearly 
point for the existence of two different regimes of growth, 
with growth exponents /3i and 02, when only two species 
of particles are deposited. The observation of the data 
for the time evolution of the roughness shows a change 
in its behavior, from an uncorrelated growth, whose ex- 
ponent is /3i to another one, where lateral correlations 
are present, which is characterized by a second growth 
exponent, /?2- 

At the very initial times, the surface is free of lateral 
correlations, once the particles are incorporated to the 
substrate, following the rules of the random deposition 
model. The exponent /3i is near 1/2, which is the typ- 
ical value for the random deposition model. However, 
at intermediate times, height-height correlations develop 
and we find a smaller value for the growth exponent /?2- 
Finally, at long times, the roughness enters into a sat- 
uration regime, and the surface is characterized by the 
roughness exponent a. 

In Fig [2] we show the log- log plot for the time evolution 
of the roughness, where simulations were performed with 
particles selected from the Poisson distribution, Eq. (3) , 
with x = 5 and n ranging from 1 to 9. The results we 
are showing, represent averages considering 10 3 different 
samples for the linear sizes L = 1024, 2048 and 400 sam- 
ples for L = 4096 and 8192. In this log-log plot we see 
two different linear regimes at the initial times, where we 
find the value j3\ — 0.503 ± 0.005 for all the lattice sizes, 
which is close to the expected value for the random de- 
position model, whose exact value is 1/2 in any spatial 
dimension. 
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FIG. 3: Plot of the growth exponent P2 versus 1/L. The 
straight line is the best fit of the P2, for several values of L and 
surface formed by adding particles selected from the modified 
Poisson distribution, Eq. (3). 



the nonlinear model, with conservative dynamics and 
non conservative noise in (1+1) dimensions, which is de- 
scribed by the Villain-Lai-Das Sarma equation [23[ . The 
exponents of this model were exactly determined from 
the renormalization group recurrence relations, where the 
nonlinear term determines the scaling behavior and the 
upper critical dimension is d c = 4. For d = 1, the expo- 
nents are: a = 1, j3 = 1/3 and z = 3. In the deposition 
problems, we say that a given relaxation process is con- 
servative if it does not change the number of particles 
in the system. On the other hand, a noise is said to be 
nonconservative when its correlation function at coordi- 
nates (x,t) and (x',t') is of the white noise type, that 
is, it is the product of two Dirac's delta functions: one 
of them, is a function of the difference between the two 
spatial positions (x — x'), and the other a function of the 
difference between the two times, (t — t'). 



At intermediate times, the interface width grows more 
slowly, with another growth exponent, ,82, which is not 
in the same universality class of the random deposition 
model. Figure [2] shows that this exponent depends on 
the lattice size L, specially for the smallest system sizes, 
which is a typical finite size behavior. For larger system 
sizes P2 smoothly increases with L, reaching a limiting 
value when L — > 00. In Figure [3] we show the plot of the 
growth exponent /?2 as a function of 1/L. The extrapo- 
lation to L — > 00 gives the best value for this exponent, 
@2 = 0.310. Fortuitously, it is close to the value 1/3, 
which is the exact result for the ballistic deposition in 
(1+1) dimensions, described by the KPZ equation. 

On the other hand, at long times, when the surface 
width reaches its saturation value for each lattice size, 
we can estimate through equation [5] the value of the cor- 
responding roughness exponent. We have found that a 
= 0.94+0.01. Another method to determine a is by us- 
ing the rms fluctuation of the squared roughness in the 
steady state, presented in Refs. [2l|,|2^. This procedure 
gives an estimate for the roughness exponent, which is 
much less dependent on the finite-size corrections. The 
effective roughness exponent is defined as 

l ln[a(L)/a(L/2)} 

a{L) = 2 ha ' (5) 

where a is defined by 



0" = \J < Wsat 2 > - < Wsat > 2 ■ (6) 

Using this method we have found that a = 0.942, which 
is in good agreement with our estimate based on the 
equation [21 With the values of a and /?2, the dynam- 
ical critical exponent is z = 3.00, which comes from the 
Family- Vicsek scaling relation. After deposition, parti- 
cles can not migrate either downwards or upwards, and 
there is no surface diffusion in the present model. The 
exponents we calculated, a, fi% and z fit very well to 



IV. IDENTICAL PARTICLES 

For surfaces formed by depositing identical particles, 
only particles with size in the range 2 to 9 lattice units of 
the linear substrate were considered. In our simulations, 
the particles are fiat, with variable horizontal length (2 
< N < 9 ) and vertical height corresponding to one unit. 
The deposition rules are exactly the same as the ones 
presented in the previous Section. The dependence of 
the surface width on the size of the particle is very weak 
and its behavior as a function of time for different lattice 
sizes is very similar to that seen for the deposition of 
particles of different sizes. We also find two different 
growth regimes characterized by the exponents fi\ and 
/?2- In Table I we display for each particle length (N), the 
values of the exponents j3% , /?2 » o, and z for the lattice size 
L = 8192. Then, even when the particles are of the same 
size, we find the this model is in the same universality 
class of the nonlinear model with conservative dynamics 
and non conservative noise in (1+1) dimensions. 



Particle size (N) 


L = 8192 




Pi 




a 


z 


2 


0.50±0.04 


0.32±0.02 


0.95±0.01 


2.98 


3 


0.50±0.02 


0.32±0.02 


0.93±0.01 


2.94 


4 


0.51±0.03 


0.32±0.01 


0.94±0.01 


2.94 


5 


0.51±0.07 


0.32±0.03 


0.94±0.01 


2.92 


6 


0.50±0.03 


0.31±0.04 


0.95±0.02 


3.05 


7 


0.50±0.01 


0.31±0.04 


0.97±0.01 


3.07 


8 


0.51±0.04 


0.32±0.03 


0.95±0.01 


2.96 


9 


0.52±0.03 


0.33±0.02 


0.95±0.01 


2.89 



TABLE I: The growth exponents for different paticles ' sizes. 



V. POROSITY 

There are few studies in the literature considering the 
formation of voids inside the volume [l6l . [22j , and only 
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FIG. 4: Porosity versus time, for a bulk formed by depositing 
particles of different sizes on a linear lattice of size L = 4096. 
The particles are chosen from the modified Poisson distribu- 
tion, Eq. (3). The inset shows the porosity at the initial 
times. 



deposits formed by the deposition of particles with two 
sizes have been analysed. In this section we present the 
results for the global porosity, where particles of differ- 
ent sizes are taken into account. The porosity, defined 
by Eq. |4]and calculated at each Monte Carlo step, sat- 
urates at the early stages of deposition, and exhibits a 
strong dependence only on the size of the particles. For 
deposits formed by a mixture of particles, selected from 
the modified Poisson distribution, or by identical parti- 
cles with size in the range 2 < N < 9, we didn't observe 
any dependence of the porosity on the size L of the linear 
substrate, for lattice sizes larger than L = 512. 

Fig. H] shows the time evolution of porosity for a de- 
posit formed by a mixture of particles, whose sizes were 
selected from the modified Poisson distribution, Eq. (3) , 
with x = 5 and n ranging from 1 to 9. The particles 
are aggregated on a linear substrate of size L = 4096, 
and the resulting deposit is a two-dimensional structure. 
Just after 50 time units we reach the stationary porosity 
value P = 0.445 ± 0.004. 

When we consider the formation of pores due to the 
deposition of identical particles with size larger than one 
lattice spacing, we find that the stationary porosity in- 
creases with the size of the particle. The deposition of 
particles with size exactly equal to one lattice spacing 
of the substrate does not lead to the formation of pores 
because the deposition follows the rules of the random 
deposition model. The inclusion of a very small fraction 
of particles of size 2 is sufficient to generate a porous 
structure, which percolates over the whole deposit [HI- 

As to be expected, the porosity increases with the size 
of the deposited particle. When a large particle is incor- 
porated to the substrate, a shadow zone is created below 
it as we have seen in Fig. [T] In this way, we expect the 
porosity increases with the size of the particle. We show 
in Fig. [5] the behavior of the stationary porosity as a 
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FIG. 5: Porosity versus particle size N in a lattice with L = 
4096. 

function of the size of the deposited particle. 

We note that the increase of the porosity is large for 
small particle sizes, but it reaches a value near its sat- 
uration for particles of intermediate size. It is easy to 
understand that the saturation value for this problem is 
P = 0.5. A simple reasoning shows that the maximum 
value of P is 0.5. This happens when the size of the par- 
ticle to be deposited is N = L/2, i.e., half of the lattice 
size. After a particle of this size is incorporated to one of 
the ends of the substrate, the probability that the next 
particle is deposited at the same height as the former 
one is about 1/ L, which is a very small number when we 
take larger substrates. On the other hand, if the particle 
is not deposited exactly at one of the ends of the linear 
substrate, no other particle can be incorporated with the 
same height. Then, in this case, only a single particle is 
deposited at each Monte Carlo step and we have P = 1/2. 
For values of L/2 < N < L, the porosity decreases, and 
is given by 



and only a single particle is aggregated in each Monte 
Carlo step, which happens, in fact, in the first trail of 
deposition. Our results show that, even for particles of 
small size, as for instance N = 9, wc find that P = 
0.489 for L = 8192, which is very close to the maximum 
possible value, P = 0.5. 

VI. CONCLUSIONS 

We studied a surface growth model where particles of 
different sizes are deposited on a linear substrate. Parti- 
cles with sizes in the range 1 < N < 9 are selected from a 
modified Poisson distribution and fall from random posi- 
tions over an initially fiat surface. We also considered the 



6 



deposition of identical particles with size larger than one 
lattice spacing of the linear substrate. Through Monte 
Carlo simulations we found that the surface width evolves 
in time following there different behaviors. At the initial 
times the surface roughness behaves as in the random 
deposition model, and the growth exponent is /3\ = 0.50. 
At intermediate times, the surface roughness grows more 
slowly with the exponent /?2 = 0.310. Finally, at long 
times, it enters into the saturation regime, character- 
ized by the roughness exponent, a = 0.94. The esti- 
mate dynamic exponent is z = 3.0. These figures put the 
present deposition model in the same universality class 
of the nonlinear model with conservative dynamics and 
non conservative noise in (1 + 1) dimensions. The two- 
dimensional bulk formed by the deposition of particles 



larger than one cell unit of the linear substrate reveals us 
a porosity increasing very fast at the initial times, and 
also reaching a saturation value, which depends on the 
size of the particle. For particles of size N = 9 the sta- 
tionary porosity reaches a value which is 2% away of the 
maximum possible value P = 0.5. 
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